Supersymmetric factorization yields exact solutions 
to the molecular Stark effect problem for "stretched" states 
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the | J = m,m) family of "stretched" states. The analytic expressions for the wavefunction and 
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In 1983, Gcndcnshtcin demonstrated that 
Schrodinger's equation is exactly solvable if the po- 
tential and its superpartner exhibit shape-invariance [T]. 
Whereas a supersymmetric Hamiltonian can be 
constructed for any potential whose ground-state wave- 
function is analytic, shape invariance only exists for 
supersymmetric potentials that are interconvertible by 
a change of a parameter other than the integration 
variable itself [3] . Herein, we make use of the methods 
of supersymmetric quantum mechanics to arrive at exact 
wavefunctions and other eigenproperties of molecules 
subject to nonresonant electric fields in closed form. 

In our previous work on the molecular Stark effect, we 
showed that for polar molecules, combined collinear elec- 
tric and nonresonant radiative fields can synergetically 
produce spatially oriented pendular states, in which the 
molecular axis librates over a limited angular range about 
the common field direction [H[S]. These directional states 
comprise hybrids of the field- free rotational states | J, to) , 
with a range of J values but a fixed value of to, which re- 
mains a good quantum number by virtue of the azimuthal 
symmetry about the fields. This has proved an effec- 
tive and versatile means to produce oriented molecules 
for applications ranging from molecule optics and spec- 
troscopy to chemistry and surface science [6HS]. How- 
ever, the eigenproperties of the Stark states in question 
could only be found numerically, typically by diagonaliz- 
ing a truncated Hamiltonian matrix. This renders all the 
treatments of molecular processes in fields, such as cold 
collisions or collective behavior of ultracold polar gases, 
analytically unsolvable |9|- Here we show that super- 
symmetric factorization of the Hamiltonian yields exact 
wavefunctions | J = to, m; lo, Alo) in closed form for a par- 



ticular ratio of the parameters lo and Alo that determine 
the interaction strengths of the molecules with the static 
and radiative fields, respectively. We found that, in semi- 
classical terms, this ratio originates in the integrability of 
the differential equation for the system's action. We also 
found that in the exactly solvable field-free and strong- 
field limits the supersymmetric problem indeed exhibits 
shape-invariance. 

We consider a 1 S molecule with a rotational constant 
B, a permanent dipole moment fj, along the intcrnuclear 
axis, and polarizability components ay and a± parallel 
and perpendicular to the internuclear axis. The molecule 
is subjected to an electrostatic field £ combined with a 
nonresonant laser field of intensity /, whose linear po- 
larization is collinear with e. With energy expressed 
in terms of B, the Hamiltonian takes the dimensionless 
form g], 

H = J 2 - lo cos - (Aw cos 2 8 + , (1) 

with the dimensionless interaction parameters lo = fJ,s/B, 
Alo = cj|| — and lo\\j_ = ^ a \\,±I/{Bc). 

The common direction of the collinear electrostatic and 
linearly polarized radiative fields defines an axis of cylin- 
drical symmetry, chosen to be the space-fixed axis Z. 
The projection, to, of the angular momentum J on Z 
is then a good quantum number while J is not. How- 
ever, one can use the value of J of the field-free ro- 
tational state, Yj im (9,(j)), that adiabatically correlates 
with the hybrid state as a label, designated by J, so 
that | J, m;Lo, Alo) — > Yj. m for to, Alo — > 0. For arbitrary 
values of the interaction parameters u> and Alo, the so- 
lution to the Schrodinger equation with Hamiltonian (fTl) 
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is an infinite coherent superposition of the field-free ro- 
tor wavef unctions, whose expansion coefficients can be 
obtained by truncating the series and diagonalizing the 
Hamiltonian in the resulting finite basis set. 

The axial symmetry of the problem allows to separate 
angular variables and express the dependence on the az- 
imuthal angle <f> via the good quantum number m. The 
Schrodinger equation for Hamiltonian then becomes 
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We note that all rotational levels are uniformly shifted 
by In what follows we use E = E^ tOL + lj±_ instead 
of the 'true' molecular energy, E^ a . Moreover, since 
the Stark effect does not depend on the sign of to, we 
define the projection of the angular momentum on Z as 
a positive quantity, to = \m\. 

By means of the substitution, if>(6) = f(9)(sia.6)~i, 
eq. ([2| can be transformed to a one- dimensional form [10] , 
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„ 4 — ui cos 8 — Aoj cos 2 8 — 
i a 6 4 



f(8)=Ef(8), 
(3) 

which will be shown to play the role of one of the requisite 
superpartner equations leading to ground state energy 
E = E . 

In what follows, we invoke supersymmetry to find an- 
alytic solutions to eq. ^ and subsequently to eq. 
Supersymmetry makes use of the first-order differential 
operators, A ± = =F^ + W(8), with W{8) being the su- 
perpotential. The superpartner Hamiltonians are defined 
by 



A ± A T = -- 



dd 2 



(4) 



q(9) can be determined, leading to the following expres- 
sion for W(8), 



W{8) 



cot 8 + (3sm8, 



(5) 



and the corresponding SUSY partner potentials, 

v!_ 1D \e) = ^-2/3(to + 1)cos(9-/3 2 cos 2 (9-£;o-t 

sin 8 4 

(6) 



rin-i , , (m + l) 2 — \ to 1 

V+ (6>) = ^ — i - 2/3to cos 8 - (3 2 cos 2 6 - E - - 



sin 2 8 



(7) 

with Eq = to(to+ 1) — f3 2 . The strengths of the combined 
fields are connected with f3 via the following expression, 



Aw = 



4(m + l) 2 



(8) 



The ground state wavefunction f (8) can be obtained 
from superpotential ^ in closed form [TTj . 

f a -{8) = 7V(-l) m (sin6l) (m+1/2) e^™ 6 ; (9) 

the normalization constant N can be expressed analyt- 
ically via the hypergeometric functions [TH [15]. The 
phase factor (—1)™ leads to the correct asymptotic be- 
havior of the wavefunction (8) (sin which, for 
j3 = 0, reduces to the ground state wavefunction of a rigid 
rotor with J = to, Y m , m (8,Q). Since the ground-state 
wavefunction ([£]) is normalizable and obeys the annihila- 
tion condition A~ = 0, the supersymmetry obtained 
is unbroken [TQ US HZ] . 

The 3D potentials corresponding to a molecule inter- 
acting with the combined fields, can be obtained by com- 
paring the K| 1D) and yi 1D) , eqs. with eqs. ^ 
and pL The detailed analysis of the ID and 3D super- 
symmetric partner potentials is presented elsewhere [14] . 

Eq. ([7]) coincides with the potential for a rigid rotor in 
the combined fields whose projection quantum number is 
to + 1 and whose interaction strengths are related by 



with the one-dimensional partner potentials V^ 10 " 1 (8) = 
W 2 (8) ± W'(9). The superpartner Hamiltonians have 
the same energy spectra except for the ground state, 
i.e., E~ +1 = E+ and Eq = 0. If the eigenfunctions 
of one of the partner Hamiltonians -H=p are known, the 
eigenfunctions of the other can be obtained analyti- 
cally via the intertwining relations, ip^-i ~ ^"V'n an d 

vv-a+vvU mug. 

For a molecule in the combined fields, we assume the 
superpotential to have the form W(8) = a cot 8 + q(9), 
where the first term corresponds to the field-free ro- 
tor 0G3 H3]. By identifying the effective potential in 
eq. (ij) with V_[ 1D \9), the constant a and the function 



Auj = 



Am 2 



= /3 2 



(10) 



Hence, given a value of (3, the Hamiltonian of a molecule 
with a projection to of the angular momentum on the 
combined fields whose interaction parameters arc related 
via eq. ^ has the same set of eigenvalues as a molecule 
with a projection to + 1 on the combined fields whose 



interaction parameters are related by eq. (10 1. 

Unfortunately, the partner potentials (J6j) and §f§ are 
not shape invariant and so eq. (|3| for a molecule in the 
combined fields is not, in general, exactly solvable [4"I|15|. 
We also note that none of the known shape-invariant su- 
perpotentials listed, e.g., in refs. [HUH] leads to exactly 



FIG. 1: Correlation diagram for a molecule in combined fields 
as a function of /3. Eigenstates that can be obtained in closed 
form are shown in red. The equidistant energy levels in the 
strong-field limit, /3 — > oo, are labeled by the 3D librator 
quantum number v. See text. 




solvable partner Hamiltonians that can be experimen- 
tally implemented for molecules in nonresonant fields. 

Figure [T] shows the energy levels of a molecule in com- 
bined fields for different values of the field-strength pa- 
rameter j3. In the weak-field limit, j3 — > 0, the energy 
levels approach those of a free-rotor, which is solvable 
exactly for all the eigenstates. For nonzero but weak 
fields, (3 = 1, the levels split into J + I components 
due to the Stark effect. In this case, the SUSY partner 
Hamiltonians are not shape- invariant, and the problem 
is analytically solvable only for the "stretched states," 
with J — m. With increasing interaction strength, the 
energies of the stretched states come closer to one an- 
other and, in the strong- field limit, j3 — > oo, coalesce into 
the ground state level of the 3D harmonic librator. In 
the strong- field limit, the supersymmetric problem be- 
comes shape-invariant again, and is exactly solvable for 
all eigenstates in closed form; the equidistant levels (la- 
beled by the quantum number v) are infinitely degenerate 
and separated by an energy difference of (2oj + 4AU;) 1 / 2 . 
The weak-field and strong-field limits are described in 
detail in a forthcoming paper |14) . 

Exact solutions for molecules in combined fields allow 
to derive molecular properties analytically. The space 
fixed dipole moment, fj,z, is given by the orientation co- 
sine, (cos 6*) = (ip(0) \ cos 9\ip(9)), and for the exact wave- 
function of eq. ([9]) can be evaluated in closed form, 

Hz/fi = (cos 9) = (11) 



FIG. 2: Ground-state energies Eq (in units of the rotational 
constant B), space fixed dipole moments Hz/h = {cos 6), 
alignment cosines (cos 2 0) , and expectation values of the an- 
gular momentum (J 2 ) for different [J = m,m;/3) states as a 
function of the interaction parameter /3. 

where I n (z) is a modified Bessel function of the first 
kind jinj . Fig. [2] shows spaced fixed dipole moments cor- 
responding to the states | J = to, to; 0) for several values 
of to as a function of the j3 parameter. The value of \iz 
rapidly increases with f3. For instance, for to = 0, it rises 
from only 0.54/1 at /? = 1 to 0.83// at /3 = 3. In the case 
of the much studied 40 K 87 Rb molecule, which possesses 
a dipole moment /1 = 0.589 Debye and a polarizability 
anisotropy Aa = 54.21 A 3 [El [20], relatively weak fields 
of e = 38 kV/cm and / = 1.75 • 10 9 W/cm 2 (correspond- 
ing to = 5) give rise to a strongly oriented ground state 
with nz = 0.9/i. This value of (cos 9) corresponds to the 
molecular axis confined to librate within ±26° about the 
common direction of the fields. 

The alignment cosine, (cos 2 9) = (ip(9)\ cos 2 9\ip(9)), 
characterizes the molecule's alignment along the Z axis 
and takes the analytic form, 

2 2/3 2 A(; m + 7/2; /3 2 ) + p^Q ™ + 5/2; /3 2 ) 

(COS u) — ~ , 

2 F 1 (;to + 3/2;/3 2 ) 

(12) 

with qF\{; a; z) — o-Fi(; a; z)/T(a) a regularized confluent 
hypergeometric function [T5]. Fig.[2]shows (cos 2 9) of the 
states I J = to, to; j3) for several values of to as a function 
of the /? parameter. The |0,0;/3) state exhibits quite a 
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strong alignment with the alignment cosine rapidly ap- 
proaching with increasing (3 the value of 0.8, which cor- 
responds to a libration of the molecular axis about the 
polarization vector of the radiative field with an angular 
amplitude of 27°. 

The expectation value of the angular momentum is re- 



lated to the orientation cosine, eq. (Ill, via 



(J 2 



.1 



m+3/2 



(2/3) 



In 



-1/2 



(2/3) 



+ 0{coaO), (13) 



We note that the dependence of (J 2 ), shown in Fig. [2] 
becomes asymptotically linear in /3 for all the values \m\, 



cf. eq. ([13). 

By making use of eqs. ^ and ([8]), one can show that 
the eigenenergy becomes: 



E,, = (J 2 ) +m 2 ( ) -2/3(m + l) (cos 0) - 1 (cos 2 0) 

(14) 

Hence by invoking supersymmetry, we found a condi- 
r , under which the 3D molecular Stark 



tion, Aw = . 
effect problem for a polar and polarizable molecule sub- 
ject to to a combination of collincar nonresonant electric 
fields becomes exactly solvable for the \ J — m,m;u>, Aui) 
family of "stretched" states. We also demonstrated 
that SUSY factorization enables constructing families of 
isospectral potentials, realizable with combined fields. 

Interestingly, it is possible to glean the reason as to 
why the exact solution of eq. ([3| is obtained for only 
one particular relation between the field strength param- 
eters, eq. p]), from the semiclassical (WKB) approxima- 
tion [21] . The eigenfunction of a ID Schrodinger equation 
assumes the WKB form f(0) oc exp[iS(0)], with S{0) the 
action of the underlying classical system. For the ground 
state of the potential the action satisfies the differ- 
ential equation, S'(0) 2 — iS"{0) = V_[ 1D \ whose solutions 
are obtained by expanding S(0) in powers of h. It turns 
out that in the case of the combined field strengths con- 
nected via eq. ([8]), the series converges to the following 
exact expression, 



S(0) 



2i 



[2/3 cos + (2m + 1) ln(sin t 



(15) 



which, when substituted into f(0), yields the exact, 
closed form wavefunction, eq. 

We note that the exact |0,0;/3) wavefunction can be 
also obtained as a "curious eigenproperty" by the method 
outlined by von Neumann and Wigner in 1929. They 
showed that by imposing the integrability condition on 
the sought wavefunction, a class of potentials could be 
derived that support a localized bound state embedded 
in the continuum [22J, [23] . 

The analytic expressions for the characteristics of the 
strongly oriented and aligned states provide a direct ac- 
cess to the values of the interaction parameters required 



for creating such states in the laboratory. Moreover, the 
available analytic eigenproperties could serve to simplify 
and, simultaneously, render more accurate, models of 
many-body systems subject to electric fields, a common 
scenario for, e.g., ultracold polar gases [2"4"] . 
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